Mon. Not. R. Astron. Soc. 000, dini (0000) 



Printed 2 February 2008 



(MN WF^ style file v2.2) 



Third Dredge-up in Low Mass Stars: Solving the LMC 
Carbon Star Mystery 

Richard J. Stancliffe*\ Robert G. Izzard^ and Christopher A. Tout^ 

^Institute of Astronomy, The Observatories, Madingley Road, Cambridge CBS OHA 
^ Carolune Institute for Quality Astronomy, http://freya.phy s.uniromal.it/~rgi/ciqua/\ 

Accepted 0000 December 00. Received 0000 December 00; in original form 0000 October 00 



ABSTRACT 

A long standing problem with asymptotic giant branch (AGB) star models has been 
their inability to produce the low-luminosity carbon stars in the Large and Small 
Magellanic Clouds. Dredge-up must begin earlier and extend deeper. We find this for 
the first time in our models of LMC metallicity. Such features are not found in our 
models of SMC metallicity. 

The fully implicit and simultaneous stellar evolution code STARS has been used to 
calculate the evolution of AGB stars with metallicities of Z = 0.008 and Z = 0.004, 
corresponding to the observed metallicities of the Large and Small Magellanic Clouds, 
respecitively. Third dredge-up occurs in stars of 1 Mq and above and carbon stars 
were found for models between 1 Mq and 3 M© . We use the detailed models as input 
physics for a population synthesis code and generate carbon star luminosity functions. 
We now find that we are able to reproduce the carbon star luminosity function of 
the LMC without any manipulation of our models. The SMC carbon star luminosity 
function still cannot be produced from our detailed models unless the minimum core 
mass for third dredge-up is reduced by 0.06 Mq. 
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• 1 INTRODUCTION 

Carbon stars are stars which show the features of carbon 
rich molecules such as C2 and CN in their spectra indicat- 
ing they have a surface carbon-to-oxygen ratio (by number) 
greater than unity. They are typically identified in photo- 
metric surveys using the colours from narrow band filters 
centred near 7800A and 8IOOA (Cook & Aaronson 1989). 
The first filter focuses on a TiO absorption feature whilst 
the latter centres on a CN absorption feature. A census of 
carbon stars has been performed in the Large Magellanic 
Cloud (LMC) and Small Magellanic Cloud (SMC). In the 
LMC a total of 7750 carbon stars have been found; in the 
SMC 1707 are known (Groenewegen 2004). The distances to 
the LMC and SMC are well determined and so we are able to 
construct luminosity functions for the carbon stars in both 
these locations. We are thus presented with an opportunity 
to test our models of stellar evolution. 

It is believed that two populations of carbon stars exist 
which reflect two different formation routes. These popula- 
tions are referred to as intrinsic and extrinsic. The intrinsic 
carbon stars are believed to be thermally pulsing asymp- 
totic giant branch (TP-AGB) stars. It has long been known 
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that the unstable double-shell burning in AGB stars leads 
to the phenomenon of third dredge-up (TDUP). The ashes 
of helium burning are mixed into the convective envelope 
enhancing the surface abundance of carbon (Iben & Ren- 
zini 1983). In low-mass stars repeated occurence of TDUP 
gives rise to carbon stars because the products of helium 
burning are carbon plus a little oxygen. However, for stars 
of around 4 M0 or greater (the limit varying with metallic- 
ity) , the temperature at the base of the convective envelope 
can become high enough for the CN cycle to operate. Car- 
bon is converted into nitrogen via proton captures (e.g. Iben 
& Renzini 1983). This prevents the conversion to a carbon 
star. If sufficient mass is stripped from the envelope of such 
a star then the hot-bottom burning can cease and a carbon 
star can then form. The extrinsic carbon stars are believed 
to be formed by the accretion of carbon-rich material from 
a more evolved donor (Van Eck et al. 1998). 

There is a long-standing problem with forming intrin- 
sic carbon stars. Detailed stellar evolution models have so 
far proved unable to produce dredge-up at low enough core 
masses and hence at low enough luminosities. When Iben 
formulated the carbon star mystery (Iben 1981) calculations 
showed that dredge-up did not occur for core masses below 
about 0.6 M0 (e.g. Sackmann 1980, Wood & Zarro 1981). 
This is also borne out by more recent calculations but leads 
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to problems in reproducing the observed luminosity func- 
tions of the LMC and the SMC. Using the results of Karakas 
et al. (2002), Izzard et al. (2004) determined that the core 
masses at dredge- up would have to be lower by 0.07 M©. 

Most modern stellar evolution codes do not solve the 
equations of stellar structure, nuclear burning and mixing 
simultaneously. Mixing is typically solved for after the struc- 
ture and burning have been determined. Examples include 
the codes used by Straniero et al. (1997) and Herwig (2000). 
Karakas et al. (2002) employ a partial simultaneous ap- 
proach. A structure iteration is followed by a mixing iter- 
ation before another structure iteration is made. In this way 
a model is converged. Initial work on a fully implicit and si- 
multaneous solution of these equations (whereby the struc- 
ture, burning and mixing are all treated together in the same 
iteration) together with an improved treatment of diffusive 
mixing, on the TP-AGB was done by Pols & Tout (2001). 
More recently StancUffe, Tout & Pols (2004) found that it 
was possible to obtain dredge-up at lower core masses than 
was previously thought. 

Following on from that work, the evolution of stars of 
between 1 and 6 M© at Z = 0.008 and Z = 0.004 is inves- 
tigated. The outcome of this detailed modelling is used as 
the input physics to a population synthesis code in order 
to compare the luminosity functions predicted by the data 
with the observed LMC and SMC carbon star luminosity 
functions (CSLF). In section 2 we discuss the stellar evolu- 
tion code used and the details of the models. The results of 
the these calculations are presented in section 3. The popu- 
lation synthesis code is briefly discussed in section 4 where 
we present our luminosity functions. 



from the pre-main sequence up to and along the TP-AGB. 
No convective overshooting or mass loss was considered at 
any stage of the evolution. We use the formalism of Bohm- 
Vitense (1958) for convection with a — 1.925, based on cali- 
bration to a solar model. Models were evolved with 199 mesh 
points up to core helium burning. At this point the resolu- 
tion was increased to 499 mesh points in order to facilitate 
the transition to 999 meshpoints just before the onset of the 
TP-AGB. 

2.1 The Core Helium Flash 

Stars of masses up to around 2.3 Mq are expected to un- 
dergo a core helium flash. However, this phase of evolution 
is numerically demanding and the STARS code is currently 
not suitable for a calculation of the evolution through it. In- 
stead a model of the desired mass is run from the pre-main 
sequence up to the helium flash. The hydrogen exhausted 
core mass and the envelope composition are recorded. A 3 
Mq model is then evolved from the pre-main sequence up 
to the point where helium ignites in the core. During this 
evolution helium burning reactions are allowed to produce 
energy but not consume helium. Once helium has ignited 
under non-degenerate conditions mass is stripped from the 
envelope and the core is allowed to grow until the model 
has the desired envelope mass and core mass. The envelope 
composition is then set to that of the pre-flash model. We 
are therefore assuming that the helium flash proceeds so 
rapidly that the core mass doesn't change and that there 
is no change in the envelope composition. These are both 
standard assumptions. 



2 EVOLUTION CALCULATIONS 

The stellar evolution code STARS developed by Eggleton 
(1971) and much updated (e.g. Pols et al. 1995) is a fully 
implicit code that solves the equations of stellar structure, 
nuclear burning and diffusive mixing simultaneously. This 
is notably different to the way most codes work with mix- 
ing usually treated separately. We use the version of the 
code developed by Stancliffe, Tout and Pols (2004) with its 
modifications for use on the TP-AGB. It has an improved 
treatment of mixing and a viscous mesh to aid numerical 
stability. 

It is well known that the phenomenon of TDUP is sen- 
sitive to the numerical treatment of the convective bound- 
ary (Frost & Lattanzio 1996). It must be ensured that the 
convective boundary remains stable (see Mowlavi 1999 for 
a detailed discussion). This is achieved in the current code 
by the use of an arithmetic mean for the diffusion coffecient 
(see Pols & Tout 2001 for full details). The adaptive mesh of 
the code provides a mass resolution of around 2 x 10^^ Mq 
at the boundary of the convective envelope. This is com- 
parable with the resolution used by Herwig (2000). An in- 
crease in the resolution does not change the amount of ma- 
terial dredged-up. The code does not treat the changes in 
energy associated with mixing material of different compo- 
sition (Wood 1981). However Frost (1997) finds that if this 
physics is not included the efficiency of third dredge-up is 
reduced. 

Models of 1, 1.5, 2, 3, 4, 5 and 6 Mq were evolved 



3 RESULTS 

We are able to evolve through a significant number of ther- 
mal pulses in all our models and third dredge-up is found 
to occur in them all, including the 1 Mq models. At both 
metallicities hot-bottom burning is found to occur in models 
of 4 Mq or greater and these stars would not become car- 
bon stars without mass loss. If these stars do become carbon 
stars once hot-bottom burning has ceased they have much 
larger core masses and so contribute to the more luminous 
end of the CSLF. Eventually each evolution sequence termi- 
nates due to unresolved numerical issues but the models are 
sufficiently evolved for our purposes where we are mostly 
concerned with the early pulses of the lowest-mass stars. 

3.1 LMC models 

For the LMC models we form carbon stars in all models of 
3 Mq and below. Details of the models are given in Tabled 
Note that the formation of carbon stars in these models 
is extremely rapid. It requires only two or three thermal 
pulses with TDUP. This is because of the low abundance 
of oxygen in these models: only a small amount of carbon 
needs to be dredged up for it to become more abundant than 
oxygen. The efficiency of dredge-up is defined through the 
paramter A = AA/dup/AMc, where AMdup is the mass of 
material dredged up and AMc is the amount by which the 
H-exhausted core grows in the preceeding interpulse. The 
maximum A value of the 1 Mq model quoted in Table Q 
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Mass 


MiTP 


Mtdup 


'^max 


Ntot 


Ntdup 


Nc 


1 


0.54803 


0.57262 


0.169 


8 


6 


2 


1.5 


0.56610 


0.57127 


0.752 


40 


4 


2 


2 


0.56419 


0.57866 


0.880 


15 


3 


3 


3 


0.63039 


0.64135 


0.992 


7 


3 


3 


4 


0.81130 


0.81676 


1.100 


9 


3 




5 


0.87062 


0.87602 


1.085 


17 


3 




6 


0.95085 


0.95643 


1.098 


11 


3 





Mass 


MiTP 


Mtdup 


-^max 


Ntot 


Ntdup 


N< 


1 


0.55929 


0.58048 


0.127 


5 


5 


1 


1.5 


0.57054 


0.58849 


0.774 


12 


4 


2 


2 


0.58133 


0.60275 


0.916 


11 


3 


2 


3 


0.69939 


0.71421 


1.048 


9 


3 


2 


4 


0.82842 


0.83166 


1.097 


8 


3 




5 


0.90026 


0.90292 


1.068 


14 


3 




6 


0.96795 


0.97430 


0.954 


7 


3 





Table 1. Details of the LMC models. Mixp and M-pouP ^"^^ the 
hydrogen exhausted core masses at the first thermal pulse and the 
first pulse for which TDUP occurs respectively. The maximum 
value of A obtained during the model run is Amax- The total 
number of thermal pulses evolved through is Ntot, Nxdup is the 
number of pulses prior to the occurence of TDUP and Nc is the 
number of pulses with TDUP required to form a carbon star. Note 
that stars of 4 Mq or greater do not form carbon stars because 
of the occurence of hot-bottom burning. 



Table 2. Details of the SMC models. Mixp and Mtdup are the 
hydrogen exhausted core masses at the first thermal pulse and the 
first pulse for which TDUP occurs respectively. The maximum 
value of A obtained during the model run is Amax- The total 
number of thermal pulses evolved through is Ntot, Ntdup is the 
number of pulses prior to the occurence of TDUP and Nc is the 
number of pulses with TDUP required to form a carbon star. 
Note that stars of 4 Mq or greater do not form carbon stars. 
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Figure 1. The evolution of surface C/O with absolute bolomctric 
magnitude, measured at the time of minimum luminosity in the 
interpulse, for models of 1, 1.5 and 2 Mq. The LMC carbon star 
luminosity function peaks at M^oi = —4.9. 



is unlikely to be a true representation of the maximum ef- 
ficiency reached because numerical instabilities prevent us 
evolving this model to a steady state. 

For models of 2 Mq or below the core masses at which 
the first episode of dredge-up occurs are lower than those 
in other models. For example, Karakas et al. (2002) find 
their 1 Mq model without mass loss does not begin dredge- 
up until its core mass is 0.657 Mq. This is 0.084 Mq more 
massive than ours. This is important because we expect the 
lower-mass stars become the lower-luminosity carbon stars. 
We also find more efficient third dredge-up in these stars. 

With the results of Karakas et al. (2002), Izzard et al. 
(2004) used a synthetic evolution code to fit the observed 
CSLF. They found that the minimum core mass for third 
dredge-up to occur would have to be lower by 0.07 Mq and 
the minimum dredge-up efficiency would have to be 0.5 to 
get a good fit. Our results are in very good agreement with 
these predicitions. We find deeper and earlier dredge-up. 
Our models are also consistent with the requirement of a 
minimum temperature at the base of the convective envelope 
of 2.5 X l(f K for third dredge-up to occur as used by Marigo, 
Girardi & Bressan (1999) in their synthetic code. 

Figure shows the evolution of the surface C/O abun- 



dance by number with absolute bolometric magnitude. The 
data points are taken at the minimum luminosity occuring 
during the interpulse as this is what determines the lowest 
luminosity carbon star that forms. The results are encour- 
aging. We find that our two lowest mass models are able to 
form carbon stars at an absolute bolometric magnitude of 
about -4.2. The LMC CSLF peaks at Mboi = -4.9. 



3.2 SMC models 

As with the LMC models, we are able to form carbon stars 
in all our models below 3 Mq . The details of the models are 
presented in Table |21 Again, carbon star formation follows 
rapidly once TDUP is established. For the same reasons as 
with the LMC model, Amax is unlikely to be truely represen- 
tative of the maximum efficiency reached. 

The results of the SMC models runs are less encouraging 
from the point of view of reproducing the SMC carbon star 
luminosity function. Our low-mass models give minimum 
core masses for third dredge-up similar to those of Karakas 
et al. (2002). Our only advantage is the occurence of TDUP 
in the 1 Mq model. Figure |5| shows the evolution of the C/O 
abundance with absolute bolometric magnitude. The lowest 
mass models become carbon stars with Mboi = —4.2. This 
magnitude is only just below the peak of the SMC CSLF. 



4 POPULATION SYNTHESIS 

In order to see if we can reproduce the LMC and SMC CSLF 
we need to create a population of stars based on our detailed 
models. We utilise the synthetic TP-AGB evolution code 
of Izzard et al. (2004). The luminosity core-mass relation 
therein is found to fit our data well and so does not need to 
be refitted. As in Izzard et al. (2004), the luminosity dip after 
each thermal pulse (see Iben & Renzini 1983) is modelled 
by a factor of the form 



/l = 1 — 0.5 X min 



l,exp(-3 — ) 



(1) 



where t is the time from the beginning of the current pulse 
and Tip is the interpulse period. In this form the luminosity 
dip found in the detailed models can be reproduced (see 
Figure (SJ. It is necessary to include this dip in order to 
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Figure 2. The evolution of surface C/O with absolute bolometric 
magnitude, measured at the time of minimum luminosity in the 
interpulse, for models of 1, 1.5 and 2 Mq. The SMC carbon star 
luminosity function peaks at Mt,oi = —4.5. 
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Figure 3. Plot to of the exponential fit (solid line) to the post- 
flash luminosity dip of the 4th pulse of the 2 Mq Z= 0.008 model 
(dashed line). The fit leads to an underestimate of the bolometric 
magnitude of about 0.1 mag. 



reproduce the low-luminosity side of the CSLFs. The core 
masses at the first thermal pulse, Mitp, minimum core mass 
for third dredge-up to occur, Mc,min, and A values used by 
the synthetic code are fitted to the above detailed models. 
Additional models of 1.25, 1.75 and 2.25 MQwere made for 
the LMC in order to better model the behaviour of the core 
mass. The fits are 



MiTp = 0.9557 + 0.31741 exp 



for the LMC and 
A/iTP = 0.56633 -I- 



(A'/* - 1.164)2 
1.1192 
0.82352 

' 1 + 0.318312'^546-Af. 



0.067603M, 



1 + 0.096788 



M, -2.6354 



(2) 



(3) 



for the SMC, where M* is the mass of the star at the begin- 
ning of the TP-AGB. Third dredge-up begins in the models 
when the core mass reaches a minimum value, Mc.min, de- 



fined by 

Mc.min = 0.9659 + 0.26309 exp 



(Af, - 1.0724)2 
1.2629 
0.73763 



for the LMC and 
A^cmin = 0.58042 + ■ 



1 + 0.33442-86280-M. 

0.065376Af. 



(4) 
(5) 



1 +0.098049*^.-2-6691 
for the SMC. 

The behaviour of A once TDUP has commenced is mod- 
elled by 



A = a[l - (1 - 0.42n)e~"^], 

a = max(0, 2.0369 + 15Z - 2.8775A/* + 1.6206A^2 

+ 0.023211Af,''), 



: 6.835 



-0.335270A//; 
4.44870 



(6) 
(7) 

1 0.002"-> -2-5 ' 

where Z is the metallicity and n is the pulse number count- 
ing the first pulse where dredge-up occured as 1. The last 
factor in Equation (4) allows A to rise to a maximum over 
several pulses and then begin to fall off as is seen in the de- 
tailed models. The low number of thermal pulses calculated 
in the 1 Mq models makes a reliable estimate for A difficult 
and so for stars of between 1 and 1.5 Mq we use the val- 
ues as computed for the 1.5 Mq case. This makes dredge-up 
slightly more efficient but as carbon stars form rapidly at 
such low masses this does not significantly affect the results. 

A population of 10,000 stars was evolved for 16 Gyr. A 
Kroupa, Tout & Gilmore (1993) initial mass function and a 
constant star formation rate were assumed. Mass loss was 
included as the Vassiliadis & Wood (1993) type prescription 
used by Karakas et al. (2002). The superwind phase is turned 
on when the Mira pulsation period of the star reaches 500 d. 
We have fitted the synthetic code to results from detailed 
models without mass loss. Because the Vassiliadis & Wood 
formalism causes significant mass loss only in the later pulses 
and our models form carbon stars rapidly, the impact of this 
approximation on our results should be limited. 

The results of the population synthesis runs are shown 
in Figures 01 and |S] The models are normalised such that 
the peak of a model matches the peak of the corresponding 
observations. The model fits the observations of the LMC 
CSLF very well except for the model a slight underabun- 
dance of carbon stars between Mboi ~ —4 and Mboi = —4.75. 
Note that we do not fit the very low-luminosity carbon stars. 
These are likely to be extrinsic, rather than intrinsic, carbon 
stars (see Izzard & Tout 2004). 

The fit to the SMC CSLF is somewhat disappointing. 
The peak in the model luminosity function is almost 1 mag- 
nitude too bright. If the minimum core mass for dredge- up is 
reduced by 0.06 Mq as in Izzard & Tout (2004) then we are 
able to reproduce the SMC CSLF (see Figure]^. Again, the 
very low luminosity tail is expected to be due to extrinsic 
carbon stars. 

It should be noted that our models consist only of a 
single metallicity and do not account for the finite size of 
the clouds. Both these effects would broaden the luminosity 
function slightly. The work of Harries, Hilditch & Howarth 
(2003) gives examples of systems with distance moduli that 
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Figure 4. The theoretical fit (solid line) to the LMC CSLF. The 

histogram is observational data taken from Groenewegen (2002). 
The CSLF is reasonably well reproduced by the theoretical mod- 
els. 
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Figure 5. Theoretical fits to the SMC CSLF. The histogram is 
observational data taken from Groenewegen (2002). The unad- 
justed theoretical model (solid line) reproduces the shape of the 
luminosity function well, but it is too bright by almost 1 magni- 
tude. If the minimum core mass required for TDUP is reduced by 
0.06 Mq (dashed line) then the CSLF is well reproduced. 



can be as much as 0.2 magnitudes from the quoted distance 
to the SMC. This is unhkely to have a great effect on the 
luminosity function, given the binning is of this order but it 
may improve the LMC model. The effect of metalhcity vari- 
ations is much more difficuh to quantify, though we note 
that varying the metaUicity by a factor of 2 (i.e. going from 
Z= 0.008 for the LMC to Z= 0.004 for the SMC) does pro- 
duce noticeably different results. 



5 CONCLUSION 

We have made evolution calculations, with a fully implicit 
and simultaneous code, of the thermally pulsing asymp- 
totic giant branch evolution of low- and intermediate-mass 
stars at LMC and SMC metaUicities. For the LMC we find 
more efficient dredge-up at lower core masses than others. 
The SMC models are comparable to those of Karakas et al. 
(2002). Using a population synthesis code calibrated to the 
detailed models we find we are able to reproduce the carbon 
star luminosity function of the LMC. The SMC CSLF can 
still only be reproduced if the minimum core mass required 
for TDUP is reduced by 0.06 Mq. 

The fact that these calculations can reproduce, for the 
first time, the LMC CSLF is of great interest. It is not 
clear why the present calculations succeed while others have 
failed, but the simultaneous solution of structure and mixing 
may be important. We are currently exploring the reasons 
for the difference between our results and those of others. 
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